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Abstract 

A review is presented of the different theoretical models proposed to approach con- 
sistently the interplay between soft an^ 
imentally at HERA for the first time. 
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^ ! sistently the interplay between soft and hard physics, that can now be studied exper- 
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&: 1 HERA Physics 
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• rH , The range of the physics scanned by HERA is very broad, from hadroproduction to 

/S ■ 

^ • jets, charm, quark fragmentation and instantons, even though in this contribution we 

will only focus on low x. 

The setting up of HERA, the first electron-proton collider, at DESY during the 
nineties of the last century made possible the experimental study of deep inelastic scat- 
tering (DIS) processes and, more generally, of quantum chromodynamics (QCD), un- 
der kinematical conditions where the interplay between soft and hard physics should 
play an important role 

The most important measurement in DIS is that of the cross-section of the process: 

ep — > eX, 
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as a function of any pair of independent Lorentz invariants built from the kinematic 
variables g, the four-momentum transfer mediated by the virtual photon, and P, the 
four-momentum of the incoming proton, e.g. {Q'^,x): 

= - = (1) 

The large center of mass energy s provided by the ep collider allows the detectors 
of the two international experimental collaborations working at HERA, HI p| and 
ZEUS the study of the kinematic regimes both at large and very small x. This 
newly accessible region of very small x is specially interesting from the theoretical 
point of view. Since 

^= n2^^ 2 > (2) 

at fixed the limit x <^ 1 corresponds to the limit of large s, s ^ (Regge limit), 
where the Regge Field Theory, which was used to describe the hadronic processes 
before QCD was accepted as the general field theory accounting for the strong inter- 
action, should be valid. Thus, HERA makes possible the study of DIS where both the 
perturbative QCD limit, ^ A^, and the Regge limit should apply, i.e. where the 
hard-perturbative and soft-nonperturbative physics should interplay and shed light 
on the fundamental question of confinement. On the other hand, very small x means 
parton densities in which the proton momentum fraction is being shared by many 
partons (mainly gluons), i.e. very high density parton densities, and thus this region 
of very small x connects the study of DIS with the physics of the heavy-ion collisions, 
so important in the next future with the operative start of the Large Hadron Collider 
(LHC) at CERN. 

Theoretically, the DIS differential cross-section can be expressed in terms of two 
independent structure functions and F2(x,(5^): 



(Pa Ana' 



dxdQ'^ xQ^ 



(3) 



where a is the electromagnetic coupling constant, y is the fraction of the electron 
energy tranferred to the proton in the proton rest frame (0 < y < 1), and Fl{x, Q"^) 
is the longitudinal structure function: 

Fl = F2-2xF^. (4) 
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The extraction of F2 from the cross-section measurement (|^) imphes an assumption 
on Fi, since there is no direct measurements in the HERA regime. Once this 
assumption is made, we can compare the experimental F2 with the theoretical one, 
that can be written in terms of the quark and antiquark densities, qi and g^, and the 
quark charges, e^.. However, it has not been possible to derive the hadronic structure 
from first principles including the interactions of quarks and gluons as given by QCD, 
and the same happens when calculating the ep cross-section. 

In fact, only three things are rigorously calculable in QCD due to the fact 
that they involve just one large scale: the process e^e~ hadrons, F2 in DIS and 
the form factor in exclusive processes [^. In DIS, the equations that include the 
perturbative effects and give the relation between the parton distribution function 
taken at two different scales are called evolution equations, and they can be obtained 
rigorously by extracting the perturbative part of a^p in the Operator Product Expan- 
sion (OPE) frame and summing up the large logarithmic perturbative corrections in 
the Renormalization Group (RG) equations @] for the regime with a large scale: 

— » 00, with — = X ~ 1 {uglnx <^ 1). (5) 

An equivalent formulation to the one based on OPE but expressed in terms of parton 
language is the one of the DGLAP equations [0. The physical meaning of this method 
is very clear since it deals with parton densities and fragmentation functions, which 
are basic quantities for describing short- distance reactions p. From (||), we see that 
the DGLAP equations can be written by neglecting the terms in In (Q^/ s) at the 
leading log approximation in ln{Q'^/Q^). This leads to consider all ladder diagrams 
in which the transverse momenta are strongly ordered: 

Q2«---4«4^^^«---g2. (6) 

This is expected to be good enough when is large but x is not too small. In this 
approximation the DGLAP evolution equation for the parton density corresponding 
to the quark of flavor i is (a similar equation can be written for the evolution of the 
gluon density, g): 



dqi{x, Q"^) as dz 



dlnQ'^ 2ti Jx z 



(7) 
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The -P (f ) functions in (|^) are the calculable splitting functions giving the probabil- 
ities for the parton branchings. 

When X is small but is not very large the DGLAP equations can no longer be 
used. It is for this new regime: 

A' < Q2 < s, (8) 

where no theoretically rigorous evolution equation is available, that a new evolution 
equation analytically solvable at fixed has been proposed, the BFKL evolution 
equation [^], in which we can, from (P), consider the leading log approximation in 
ln{l/x), which translates in resumming the ladder diagrams with strongly ordered 

xq <ti ■ ■ ■ Xi <ti Xi+i <^ ■ ■ ■ x, (9) 

and no ordering on kxi- The BFKL equation is an evolution equation in x of the 
gluon density g{x, k^), dominant at small x: 



g{x, k'^) - g{x, k"^) ^ g{x, k^ 



k'r - ^tI Mk'S - kl 



(10) 



dln{l/x) TT Jo k'^ 

Thus g{x, k^) can be calculated for any small x once it is known at some starting 
value xq, for all k^. 

As already mentioned, the BFKL equation can be solved analytically for fixed as, 
with the X and kr behaviors of the solution being: 

^(x,4)ocf— ) 4^T-^, (11) 
\XqJ ^ Inl/x 

with 

X^o ^ (^^a,/7r)4 ln2 ~ 0.5 (12) 

for ric = 3 colors, being the famous BFKL Pomeron intercept. At this point one has 
to note that different problems appear when trying to solve the BFKL equation at 
the next-to-leading order (NLO) [l^, among them the impossibility of treating all 
the terms in the new solution as Regge terms. Still, the NLO estimates of the BFKL 
Pomeron intercept seem to lead to a value: 

A^^o ~ 0.25^0.3. (13) 
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Also attempts have been made to include the running of as(Q^) in the calculation. 

In the region where x is small and is not too large, now accessible at HERA, 
the old Regge Field Theory should work. The Regge Field Theory [ll|] is a method 



based on analyticity, crossing symmetries between the s and t channels and unitarity 
in the complex angular momenta plane, valid to compute the cross-section of hadron- 
hadron collisions at high energies (the equivalent in DIS to small x), by squaring 
the sum over the scattering amplitudes due to the mesons that can be exchanged 
in the process. These mesons can be empirically grouped along the so called Regge 
trajectories, parametrized as straight lines which relate the mass m and the spin J of 
the exchanged mesons, J = a(m^): 

a{t) =ao + a't, (14) 

where ao is the intercept and a' the slope of the trajectory. Thus the elastic cross- 
section for a hadron-hadron collision will be calculated by summing over all the Regge 
trajectories whose resonances can be exchanged in the reaction, and the final result 
can be written as: 

^ oc (/5(t))^.2"W-^ (15) 

with /3{t) an unknown real function. Now, by using the optical theorem which relates 
the total cross-section to the forward elastic amplitude, we can predict the behavior 
of the total hadron-hadron cross- sect ion: 

atot oc (16) 

The huge amount of experimental data on hadronic cross-sections for many differ- 
ent processes shows an universal and steady rise at large energies that can only be 
accounted for by parametrizing all these cross-sections as the sum of two different 
components: 

atot = + (17) 

with A and B process-dependent constants and the intercepts q;_r(0) ~ 0.5 and 
ap{0) ~ 1.08, universal process-independent constants. While the first term in 
Eq. ( p!7|) represents the exchange of the experimentally detected mesonic resonances 
(secondary Reggeon), the second term represents the exchange of a hypothetical ob- 
ject, the soft Pomeron, which has the vacuum quantum numbers (electrically and 
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color neutral, isospin and C-parity +1), and it is the one responsible, through the 
value of its intercept larger than 1 (supercritical Pomeron), for the rise of o"*°* at large 
energies. This soft Pomeron, as it also happens for the BFKL Pomeron, is suspected 
to be of gluonic nature. Of course. Equation (|17D violates the unitarity constraints 
imposed by the Froissart bound on the cross-section behavior, so we must interpret 
the soft Pomeron just as an useful phenomenological tool. 

As a matter of fact, DIS processes at small x can be viewed in terms of the Regge 
Field Theory as virtual photon-proton scattering at high energy, with F2 dominated 
by the gluon content of the proton: 



where this should be compared with the solution of the BFKL evolution equation, so 
we can write: 

F2 oc , (19) 

and we view the exchange of a hard BFKL Pomeron in the Regge language as the sum 
over graphs with one gluon ladder between the interaction particles in perturbative 
QCD. 

Now we are in conditions to address the problems at the origin of the theoretical 
models that we pretend to review in this contribution. The first question is how to 
connect QCD and the Regge Field Theory in a theory which could compute quanti- 
tatively any DIS or hadronic process in both the soft and hard regimes. The second 
problem concerns the perturbative QCD prediction of a strong increase of the parton 
densities at low x, which led to the idea |T2[ that the density of partons (gluons) 
becomes so big that at some point these gluons cannot be considered as indepen- 
dent partons any more, but they interact among them. To determine whether this 
is actually the case, and whether these interactions must be taken into account for a 
consistent description of DIS processes in this small x regime, is what is called the 
saturation problem, at present intensively analyzed, both theoretically and experi- 
mentally at HERA. To precisely establish |jl3| the way saturation relates to the idea 
of unitarity, a better understanding of the fundamental question of confinement is 
needed. 
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2 Theoretical Approaches 



Many models have been proposed in the attempt to make the connection between soft 
and hard physics in a theoretically consistent way during the last years. Given the 
current lack of both theoretical and experimental tools needed in order to articulate 
an universal and fundamental approach, all these models are phenomenological to a 
larger or smaller degree, and they can basically be classified in three main goups. The 
first class of models explicitly includes saturation, as the dipole models which use the 
idea that the virtual photon actually splits into a q-q dipole and it is this dipole which 
subsequently interacts with the proton. Here, the size of the dipole will determine 
the soft or hard character of the interaction. The second class of models, that we will 
call phenomenological parametrizations of structure functions at low Q^, provides a 
parametrization of the structure function F2{x^ Q'^) describing the experimental data 
in the (nonperturbative) region where perturbative QCD is supposed to fail, and then 
uses this parametrization to be plugged in as an initial condition in the perturbative 
evolution equations to obtain a description of the experiment in the whole kinematic 
range. The philosophy behind the third group of models (dynamical models for the 
low behavior of F2) is just the opposite. Now a QCD-based consistent description 
of the data in the perturbative region is considered, and then extrapolated down to 
the region of not large Q^, through the evolution equations and under certain ad- hoc 
conditions imposed to maintain the consistency of the description in the region where 
the perturbative treatment is supposed not to work any more. 

The first explicit realization of the idea of saturation was presented not for the 
nucleon, but for the nuclear case [Q, when it was shown that the Glauber model 



without gluon interaction leads to the saturation of the growth of the quark and 
gluon distributions at fixed s as becomes small. 



One dipole model very transparent in its physical interpretation is the GW model []T5 
Here the structure function F2 is separated into two terms: 

F2 = F^ + F^, (20) 

J d^fdz\^'''''{f,z,Q^)\^a{x,f). (21) 
7 



with, for X ^ 1: 



iT,L 



47r2 



a 



Here is the known wave function for a transverse (T) or longitudinal (L) polarized 
virtual photon 7* to split into a q-q dipole, a is the dipole cross-section describing the 
interaction of the dipole with the proton, and r is the transverse separation of the q-q 
pair. Unitarity is built in by the phenomenological form of the dipole cross-section: 

a = (To{l - exp{-r'^/4:Rl{x))) 
Roix) = il/Qo)ix/xo)^/', 
where Rq is the so-called saturation radius, Qo = 1. GeV, and parameters cxo, Xq and 
A are fitted to all inclusive DIS data with x < 0.01. Thus in this model, at small r 
one has color transparency and a strong growth of a with x: 

a ~ r^a;-^, (23) 

while at large r (or a; — 0) o" approaches the black-disk constant value ctq (saturation). 
The transition to saturation is governed by Rq{x). 

There are models including other saturation mechanisms. One of them is the in- 
teraction between partons in the parton cascade , not taken into account in QCD 
evolution equations, but that could become important to slow down the growth of 
parton densities. Here, the parton interactions will create an equilibrium-like sys- 
tem of partons with a definite value for the average transverse momentum, Qs{x) 
(saturation scale). A different proposed saturation mechanism is the percolation of 
strings [|l^], a second order phase transition which takes place when clusters of over- 
lapping strings, with size of the order of the total transverse area available, appear. 
This phase transition is used as an indication for the onset of saturation of the density 
of partons. In this approach, which can be generalized to the nuclear multiple 
exchange model for ep collisions is needed. 

For the nuclear case, the first approach ||I4[ presenting the k"^ saturation curves 
for the quark and gluon distributions at fixed s has been extended to the case 
where interaction among gluons is taken into account, obtaining the x dependence 
of those curves. Also all multiple Pomeron LO exchanges have been included [|l9l in 
deriving a small x evolution equation of F2 for a large nucleus from the first nuclear 
approach without interaction [I^l, and it is shown that in the double leading log limit 
this equation reduces to the GLR equation |12]. In a different model |2^, the valence 
quarks of the nucleons of the nucleus are treated as the sources of the small x gluon 
distribution of the nucleus. 
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Among the phenomeno logical parametrizations of structure functions at low Q^, 
the DL model is a parametrization of F2 that uses two separate (soft and hard) 
Pomerons: 

+ h{q^)x-^^{l-xy (24) 

with eo = 0.4372 (hard), ei = 0.0808 (soft), and 62 = -0.4525 (valence). The ABY 
model 1^ is another parametrization of F2 using two different components, a hard 



Pomeron plus a soft Pomeron that in this model is taken as a flat term. 

The CKMT model p3| is also a phenomenological parametrization of F2 at low Q^, 



but it uses only one effective Pomeron. The CKMT model proposes for the nucleon 
structure functions: 

F2(x, Q2) = Fs{x, Q2) + Fi,s{x, Q'), (25) 

the following parametrization of its two terms in the region of small and moderate 
Q^. For the singlet term, corresponding to the Pomeron contribution: 



1+A(Q2) 

(26) 



where the x—>-0 behavior is determined by an effective intercept of the Pomeron, A, 
which takes into account Pomeron cuts and, therefore (and this is one of the main 
points of the model), it depends on Q^. This dependence is parametrized as : 

Thus, for low values of (large cuts), A is close to the effective value found from 
analysis of hadronic total cross-sections (A~0.08), while for high values of (small 
cuts), A takes the bare Pomeron value, A~ 0.2 ^ 0.25. The parametrization for the 
non-singlet term, which corresponds to the secondary Reggeon (f, A2) contribution, 

F^six, Q') = xi--(l - x)"(«^) {-^Tb) " ' ^^^^ 

where the x— >-0 behavior is determined by the secondary Reggeon intercept or, which 
is in the range api ~ 0.4 0.5. The valence quark contribution can be separated 
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into the contribution of the u (B^^) and d (B^^) valence quarks, the normahzation 
condition for valence quarks fixes both contributions at one given value of [Ql = 
2. GeV^ has been used in the calculations). For both the singlet and the non-singlet 
terms, the behavior when is controlled by niQ"^), with n{Q^) being 

Therefore, for = 0. the behavior of the valence quark distributions is given by 
Regge intercepts, n(0)=Q;ij(0)— (0)~ 3/2, while the behavior of n{Q'^) for large 
is taken to coincide with dimensional counting rules. The total cross-section for real 
((5^=0) photons can be obtained from the structure function F2 using the following 
relation: 

(30) 

Q2=0 

The proper F2{x,Q'^)^Q'^ behavior when Q^^O, is given in the model by the last 
factors in Equations p6| ) and (^), leading to the following form of the (T^°p{v) in the 
CKMT model: 

cr;j(i/) = A-K^a [A^ a^i-^°(2mz/)^» + {B^^^ + B^^) r"«(2mz/)"«'^) . (31) 

The parameters in the model were determined from a joint fit of the 0"*°* data and 
NMC data on the proton structure function in the region 1. GeV^ < < 5. GeV^ , 
obtaining a very good description of the available experimental data. The next step 
in this approach is to introduce the QCD evolution in the partonic distributions of the 
CKMT model and thus to determine the structure functions at higher values of Q^. 
For this, the evolution equation in two loops in the MS scheme with A = 200. MeV 
was used. The results obtained by taking into account the QCD evolution in this way 
are in a very good agreement with the experimental data on ^2(0;, Q^) at high values 
of Q\ 

The ALLM parametrization ||2j 
does no introduce QCD evolution. 

Among the dynamical models of the low behavior of -F2, we should mention the 
GRV model [^], in which the QCD evolution equations are extended down to the very 
low region (Q^ < 1. GeV'^) using dynamical parton densities generated radiatively 
from valence-like inputs at some resolution scale. In this model both LO and NLO 



of F2 also uses Q -dependent powers of x, but it 
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approximations are used. Another dynamical model is the KP model [Q, where the 
DGLAP evolution equation at NLO is solved by giving analytical parametrizations 
for the parton ditributions. This model takes into account the contributions of higher- 
twist (renormalon-type) operators of the Wilson OPE, which are important at low 



. Finally, the BK model considers the contributions both from the parton 
model with QCD corrections extended to the low region and from the low mass 
vector mesons. There are other vector meson dominance (VMD) models 

On top of these, one can find models which are in the middle of two of the groups 
above. Thus, in the NZZ model a color dipole approach is used to solve the BFKL 
equation with running coupling constant, while another dipole model with both a 
hard and a soft Pomeron has been proposed pO| where the large dipoles couple to the 
soft Pomeron and small dipoles couple to the hard Pomeron. This model has been 
applied to the case of the charm structure function Fl^i^x, Q^). Also an explicit dipole 
model with the CKMT pattern of energy behavior (effective Q^-dependent Pomeron 
intercept) has recently been presented |31]. Other attempts have been made by 



interpolating between Regge behavior and the high DGLAP asymptotics [ p2| . 



3 Discussion and Conclusions 

In this contribution we have presented a basic introduction to HERA physics, mainly 
focused on the field of low x physics. The low x HERA physics provides an important 
experimental tool to obtain a better understanding of the interplay between soft and 
hard physics and the relation between concepts as saturation and unitarity, and to 
address more consistently the fundamental question of confinement. Models including 
saturation present a decrease of the effective Pomeron intercept as x ^ 0, while in 
other models the effective intercept should increase as energy (or increases. 

The experimental evidence of saturation at HERA, specifically the presence of a 
change in the dependence of F2 at very small x and moderate Q^, presented in 
a very graphic way as the turnovers of the logarithmic slopes of F2, in particular of 
dF2/dlnQ'^ at moderate and small x, is not conclusive. 

Essential information on the behavior of the structure functions in the region of 
extremely small x, not accessible at HERA, will be available at LHC in a hopefully 
near future. 
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